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ABSTRACT 

Experiments were conducted to study the effects of spray 
trajectory angles upon heat flux for flat and enhanced surface 
spray cooling. The surface enhancement consisted of straight 
fins machined on the top surface of a copper heater block. 
Spray cooling curves were obtained with the straight fin surface 
aligned both parallel (axial) and perpendicular (transverse) to 
the spray axis. Measurements were also obtained on a flat 
surface heater block for comparison purposes. Each copper 
block had a cross-sectional area of 2.0 cm 2 . A 2x2 nozzle array 
was used with PF-5060 as the working fluid. Thermal 
performance data was obtained under nominally degassed 
(chamber pressure of 41.4 kPa) conditions. Results show that 
the maximum CHF in all cases was attained for a trajectory 
angle of 30° from the surface normal. Furthermore, trajectory 
angles applied to straight finned surfaces can have a critical 
heat flux (CHF) enhancement as much as 75% (heat flux value 
of 140 W/cm 2 ) relative to the vertical spray orientation for the 
analogous flat surface case under nominally degassed 
conditions. 

Keywords: Enhancement, Spray Cooling, Finned Surfaces, 

Heat Transfer, Trajectory Spray 

INTRODUCTION 

Numerous research efforts have been undertaken to gain a 
better understanding of the general phenomena and critical 
parameters associated with the spray cooling heat transfer 
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process. Previous studies have parametrically examined the 
effect of secondary gas atomizers vs. pressure atomizers 
[17,19], mass flux of ejected fluid [4,20], spray velocity [2,16], 
surface impact velocity [2,5,15], surface roughness 
[1,11,16,17], ejected fluid temperature, chamber environmental 
conditions, and spray footprint optimization on the effective 
heat flux across the surface [10]. Other topics researched to 
date include the effect of surfactant addition [12,13], and 
secondary nucleation [9,14,17]. 

This work is a continuation of the enhanced surface study 
by Silk et al. [18], with emphasis on straight fins as the featured 
surface structure geometry. The objective of the current work is 
to examine the effects of spray trajectory angle upon heat flux 
for a flat and enhanced surface (specifically straight fins) when 
using a multi-nozzle array. 

Previous studies dealing with spray trajectory angles have 
emphasized the spray cone and its relation to the heater surface 
from either a vertical or horizontal position. The work by 
Mudawar and Estes [10] examined heat flux as a function of 
cone angle and nozzle height for a given flow rate. The heater 
surface tested was square (12.7x12.7 mm 2 ) while the nozzle 
used had a circular spray footprint. Working fluids used were 
FC-72 and FC-87. The authors determined that CHF was a 
function of volumetric flow distribution on the heater surface. 
The optimum CHF was attained by inscribing the heater surface 
with the perimeter of the spray cone footprint. 

In the work by Kearns et al. [7], spray cooling of a row of 
heaters was performed inside a narrow channel. A total of nine 
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heaters was used each with an area of 38.1 mm 2 . The channel 
had length, width and height dimensions of 400 mm, 265 mm, 
and 255 mm respectively. Fluid was sprayed into the channel at 
one end by a single full cone nozzle with a cone angle of 55°. 
The configuration was designed to simulate confinement 
conditions inside an actual circuit board. The working fluid 
used was PF-5060 at 101 kPa. The maximum dissipation 
(60 W) occurred with the lead heater closest to the nozzle. The 
minimum (20 W) occurred with the heater farthest away. The 
authors concluded that this was due to the proximity of the 
leading heater relative to the nozzle as well as the impingement 
angle on its leading edge. 

Most studies that have examined enhanced surfaces have 
done so primarily from the perspective of surface roughness. 
Sehmbey et al. [17] gives an overview of spray cooling and 
provides a comparison of its effectiveness when using liquid 
and secondary gas atomizers (air used as the secondary gas). 
Heat flux was measured and presented for both techniques. 
Both the heat flux and the convection coefficient were found to 
have comparable values for both atomizer types. The authors 
concluded that the most important parameters affecting heat 
transfer are the fluid properties, spray velocity and surface 
conditions. It was also found that the heat transfer coefficient 
increased with the use of smooth surfaces (R a < 0.1 /im) for gas 
atomized sprays, while the opposite trend was observed for 
liquid atomized sprays. 

Pais et al. [11] studied the effects of surface roughness on 
heat transfer when using spray cooling. The surface roughness 
values studied were 22, 14 and 0.3 /uel The heat exchange 
surface was copper with a projected area of 1 cm 2 . An air-assist 
atomizing nozzle was used with deionized water as the working 
fluid. Tests were conducted at a nozzle height of 23 mm. Tests 
were run to the CHF value for all surface tolerances. It was 
found that the 0.3 /mi rough surface achieved the highest heat 
fluxes of the group with a peak heat flux of 1250 W/cm 2 . 
Furthermore, the onset of nucleate boiling was experienced at 
lower superheat values for this roughness. The authors 
attributed the heat transfer enhancement to early bubble 
departure from the surface and nucleate boiling. The authors 
also concluded that secondary nucleation has a primary role as a 
heat transfer mechanism only if the surface finish is smooth. 

Studies utilizing straight finned enhanced surfaces have 
been performed in both pool boiling and spray cooling. Chien 
and Webb [3] investigated the effects of structured tunnel 
dimensions on nucleate boiling convection coefficients for heat 
fluxes ranging between 2 and 70 kW/m 2 . Tests were performed 
on a 19. 1 mm diameter horizontal tube using R-l 1 and R-123 as 
working fluids. Tunnel pitch, height, width, and base radius 
were the primary dimensions reviewed. The authors found that 
fins smaller than the height of 0.9 mm experienced significant 
increases in the convection coefficient as the fin count increased 
from 1378 fins/m to 1575 fins/m. They also found that using 
straight fins helped to increase evaporation by retaining more 
liquid. Increasing fin height had little effect upon the convection 
coefficient. Fin pitch was also seen as having little effect. 


The enhanced surface spray cooling investigation by Silk 
et. al. [18] consisted of a surface structure survey entailing 
pyramids, cubic pin fins and straight fins of similar feature 
geometry size and pitch. Tests were conducted using PF-5060 at 
nominally degassed (41.4 kPa) and gassy (101 kPa using N 2 
gas) conditions. Each of the surface structures tested showed an 
increase in heat flux relative to the analogous flat surface case 
for both degassed and gassy conditions. The straight fins 
attained the maximum values for all the surfaces tested. The 
degassed case showed a straight fin heat flux enhancement of 
55% (nominal value of 126 W/cm 2 ) and the gassy increased by 
48% (nominal value of 156 W/cm 2 ). The authors attributed 
these increases to fluid surface management. 

The present work investigates heat transfer enhancement at 
various spray angles. This study may result in techniques 
leading to higher heat fluxes at comparable thermodynamic 
conditions and spray characteristics. Droplet velocity has been 
shown as being the dominant spray characteristic impacting heat 
flux [2] for sprays normal to the heater surface. An angular 
trajectory would reduce the vertical component magnitude with 
respect to a trajectory angle normal to the heater surface, and 
may result in a reduction in heat flux corresponding to the 
velocity reduction. The current study shows that the heat flux 
can increase with angular spray trajectories. 

NOMENCLATURE 

A Area 

H structure height 

L s distance between successive structures 

P pressure 

R a Surface Roughness 

T temperature 

TC thermocouple 

V droplet velocity 

X structure feature dimension 

We Weber number 

d diameter 

e error 

h convection coefficient 

k conductivity 

/ nozzle height above heater surface 

p pitch 

q n heat flux per unit area 

u uncertainty 

jc thermocouple distance 

T weighted volume flux for concentric ring 

d spray trajectory angle 

7 straight fin planar orientation angle 

fi viscosity 

0 surface tension 

SUBSCRIPTS 

1 concentric ring 

k conductivity 
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/ 

liquid 

max 

maximum 

q 

heat flux 

surf 

surface 

T 

temperature 

V 

vapor 

X 

thermocouple distance 

32 

Sauter Mean Diameter 

0 

spray nozzle orifice 


TEST SETUP AND PROCEDURE 

The experiments were conducted using a closed fluid loop 
system. The test rig (schematic shown in Fig. 1) consists of an 
environmental test chamber, liquid pump, flow meter, micro- 
filter and a condenser. Chamber temperature and pressure were 
measured via a T-type thermocouple and a pressure sensor. 
Temperature and pressure sensors were also placed in the liquid 
line upstream of the nozzle for fluid and supply line temperature 
and pressure measurement. 

Each of the test heaters were made of oxygen free copper 
with a uniform undercoat of 2.54 /xm nickel and 1.27 /xm top- 
surface coat of gold. Heat was supplied to the test article using a 
500 W cartridge heater. The test article was placed within the 
interior of the chamber, but was separated from the excess 



Fig. 1 Spray Cooling Test Rig Configuration 

liquid by an enclosure consisting of polycarbonate housing and 
an alumina bisque ceramic top flange. The upper section of the 
copper blocks was epoxied to the ceramic flange. Temperature 
measurements in the copper blocks were sampled via five T- 
type thermocouples mounted in the upper neck of each block 
(shown in Fig. 2). Assuming steady state 1-D conduction 


through the upper neck of the block, the heat flux was 
calculated using Fourier’s Law of conduction. The surface 
temperature was determined via linear extrapolation using TCI 
and TC2. 



Fig. 2 Copper block schematic with TC locations (not to scale) 



Fig. 3 Spray Uniformity Test Schematic (not to scale) 
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and used in the positioning of the nozzle relative to the heater 
surface for repeatability. 


Prior to each test, the spray chamber and fluid loop were 
charged with PF-5060. For the degassed case, a vacuum was 
repeatedly applied to the chamber until a pressure of 41.4 kPa 
(470 ppm gas concentration) was reached. The chamber was 
allowed to attain equilibrium prior to conducting the tests. 

All tests were run at constant chamber pressure, liquid flow 
rate (200 ml/min) and constant nozzle height above the heater 
surface. Heat was supplied to the cartridge heater in increments 
of 10 W using a programmable power supply. Steady state was 
achieved at each power level before application of the next 
successive heat load. Upon dry-out (detected by a rapid increase 
in surface temperature and a rapid decrease in heat flux), power 
to the cartridge heater was shut-off. 

A Parker Hannifin prototype spray nozzle was used for each 
of the tests. The nozzle consists of a 2x2 spray cone array. Prior 
to heat flux testing, the spray nozzle uniformity was measured 
using several stainless steel tubes of varying inner diameters, a 
graduated cylinder, and a stopwatch. The largest tube had an 
inner diameter approximately the same diameter as the heated 
surface. Size, local volume flux between concentric cylinders, 
and the local volume flux between concentric cylinders 
normalized by the average volume flux over the entire heater 
surface (T) is shown in Fig. 3. A T value of unity indicates that 
the local volume flux is identical to the total volume flux 
averaged across the entire heater surface. The outer ring (A4) is 
shown as captures 60% less volume flux than the average area 
value. The volume flux is shown to gradually increase towards 
the center of the heater region. The center ring (A } ) which 
encompasses only 17.5% of the heater area in question has 
twice as much volume flux as the total for the entire area. Given 
the volume flux variation throughout the concentric rings, the 
spray may be effectively considered a non-uniform center 
biased spray for the current nozzle height and heater area. Spray 
characteristics were not investigated during this study. 

Three trajectory angles (0) were tested in this study. These 
are 0° (vertical), 30°, and 45°. The angles were taken using the 
centerline axis of the spray nozzle manifold relative to the 0° 
orientation (see Fig. 4). The nozzle manifold height (/) was held 
constant for each of the tests. Delrin spacers were fabricated 


Nozzle 

Height 


0° 30° 45° 



Fig. 4 Centerline axis trajectory survey angles 
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Cross Sectional View for Straight Fins 


Fig. 5 Enhanced Surface Geometry 


Surface 

P 

(mm) 

X 

(mm) 

Ls 

(mm) 

H 

(mm) 

Flat Surface (If) 

0 

0 

0 

0 

Straight Fins (Is) 

2 

1 

1 

1 


Table 1 Enhanced Surface Geometry Summary 



Diagonal Side 

Fig. 6 Straight Fin CCD Images 


SPRAY SPRAY 



Transverse Orientation Axial Orientation 

SIDE VIEW 



Transverse Orientation Axial Orientation 

TOP VIEW 

Fig. 7 Straight Fin Orientations 
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Inherent to this study is the surface structure of the enhanced 
heater block used for tests. The dimensions of the straight fin 
surface (Is) are shown on Fig. 5 and detailed in Table 1. 
Photographs of the enhanced surface are shown in Fig. 6 . When 
spraying onto the straight fin surface at an angle other than the 
0 ° position, the structures’ planar angle relative to the spray 
rotational axis ( 7 ) presents two cases (and ultimately bounds) 
for the present study. The two cases are the transverse (y=90°) 
and axial orientations (7=0°) as shown in Fig. 7. Studies were 
performed with the straight fin surface at both of the 
orientations. 

MEASUREMENT UNCERTAINTY 

The primary quantity of interest for these experiments is 
the heat flux. The heat flux calculation has three sources of 
error. These are the conductivity, the thermocouple locations, 
and the error in the temperature measured. The conductivity 
value used was 389 W/m K with 1% error. The data acquisition 
unit used for thermocouple measurements had a signal to 
temperature conversion accuracy of ±0.1 °C. The error in the 
thermocouple location was determined to be ±0.56 nrn 
Equation 1 was used to calculate the error for the heat flux 
values reported. 



The uncertainty in the heat flux was determined to be 3.4% at 
80 W/cm 2 which corresponds to the smallest CHF for all the 
cases tested. Pressure values recorded had an uncertainty of ±3 
kPa. Flow meter measurements were attributed an error of ±1 
ml/min. 

RESULTS AND DISCUSSION 

Heat flux as a function of the surface temperature and spray 
trajectory angle for the flat surface (If) is shown in Fig. 8 . Heat 
flux as a function of the surface temperature and spray 
trajectory angle for the straight fin surface in the transverse 
(7=90°) and axial orientation ( 7 = 0 °) are respectively shown in 
Figs. 9 and 10. The heat flux is based on the projected area of 

2.0 cm 2 , as opposed to the actual total surface area exposed to 
the fluid. 

Fig. 8 shows that the heat flux increases as the trajectory 
angle is varied from the 0=0° position. The maximum CHF of 
96 W/cm 2 (20% enhancement relative to 0 : =O o ) occurred for the 
30° case while the 45° case attained CHF at a slightly lower 
value of 92 W/cm 2 . The maximum surface temperature reached 
for both non-vertical cases was approximately 67° C. The 0=45° 
case showed very good agreement with the 0 = 0 ° case in the 
single phase convection regime (T surf < 55 °C) whereas the 
0=30° case held slightly higher heat fluxes. Multiphase effects 
become more pronounced in the intermediate regime (55 °C < 


T surf ^64 °C) with the 45° case increasing to equivalent heat 
flux performance relative to the 30° case. The curves then 
separate in the high heat flux regime (T surf >64 °C) as CHF is 
approached. 
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Fig. 8 Heat Flux as a function of surface temperature and 
Trajectory angle for flat surface 
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Fig. 9 Heat Flux as a function of surface temperature and 
Trajectory angle for Transverse straight fins 

In the straight fin 7=90° study (shown in Fig. 9), the 
maximum CHF (140 W/cm 2 ) also occurred for 0=30°. This gave 
a heat flux enhancement of 11 % relative to the 0 = 0 ° straight fin 
case and 75% relative to the 0=0° flat surface case. The 45° case 
had a CHF value (135 W/cm 2 ) slightly lower than the 30° case. 
The heat fluxes for the 0=0° and 30° straight fin cases agreed 
throughout the surface temperature domain until entering the 
high heat flux regime (T surf > 60 °C). The 0=45° case had 
slightly lower heat flux values than the 0 ° and 30° cases 
throughout the entire surface temperature domain until 



Voi. Flux = 0.016 nrVm 2 s 
Nozzle AP= 379 kPa 
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transitioning to the high heat flux regime. In the approach to 
CHF, the heat flux performance increased beyond that of the 
0=0° case. 

The straight fin y=0° study (shown in Fig. 10) also showed 
the maximum CHF (132 W/cm 2 ) occurring for the 0=30° case. 
The heat flux enhancement relative to the 0=0° straight fin and 
flat surface cases was 5% and 65% respectively. CHF for the 


160.0 



20.0 30.0 


40.0 50.0 60.0 

TsurffC] 


70.0 80.0 


Fig. 10 Heat Flux as a function of surface temperature and 
Trajectory angle for Axial straight fins 


maximum heat flux for the y=90°, 0=30° case. 

The study by Chen et. al. [2] showed that the droplet 
velocity (as opposed to the spray density or droplet diameter) 
had the largest impact upon heat transfer. Due to the density of 
the spray at the flow rate tested, spray characterization 
measurements detailing droplet size and velocity were not 
performed. However, assuming negligible velocity losses during 
travel to the heater surface, the droplet dispersion velocity at the 
orifice (V 0 ) may be considered an approximation of the droplet 
inpact velocity. This velocity can be determined using the 
relation developed by Ghodbane and Holman [5]. The final 
form for the velocity relation is shown in equation 2. V mbe is the 
upstream velocity of the fluid in the tubing. The pressure drop 
across the nozzle is AP. 


V = 


V + 

V tube ^ 


2A P_l2cr 

Pi Pdy 


11/ 2 


( 2 ) 


The Sauter Mean Diameter (SMD) was calculated using 
equation 3 [4,8]. 


, ( 

— = 3.07 

d o l 


d'„‘ \ 

t 


- 0.259 


(3) 


0=45° case was 126 W/cm 2 . The heat flux curves for the 0=0°, 
30°, and 45° cases agreed well until just prior to CHF. The 
maximum surface temperature attained for the 45° case (T^ 
=67 °C) was slightly lower than that of the 0° and 30° cases. 
Table 2 gives a summary of CHF, maximum surface 
temperature, and CHF percent enhancement relative to the flat 
surface for each of the cases. The Table 2 data is arranged in 
descending order with respect to the heat flux. The straight fins 
outperformed the flat surface for each y and 0 orientation. The 
7=90° cases for the straight fins had higher CHF values than the 
7=0° cases. Furthermore, each set of experiments showed that 
the 0-30° case attained the highest heat flux. Table 3 shows the. 
The maximum enhancement of 75% corresponded with the 


Surface 

7 

e 

q"max 

(W/cm 2 ) 

CHF 

Enhancement 
(relative to 
Flat Surface 
CHF at 0°) 

Tmax 

(°C) 

Is 

90° 

30° 

140 

75% 

70.5 

Is 

90° 

45° 

135 

69% 

70.3 

Is 

0° 

30° 

132 

65% 

69.8 

Is 

0° 

45° 

126 

58% 

66.7 

Is 

— 

0° 

126 

58% 

69.1 

If 

— 

30° 

96 

20% 

67.4 

If 

— 

45° 

92 

15% 

67.8 

If 

--- 

0° 

80 

0% H 

70.0 


Table 2 Summary of Trajectory Study Test Data 


The We (shown in equation 4) also is of particular 
importance. The We (which incorporates a second order 
velocity term) is the ratio of the inertia to surface tension forces. 
It is a metric for the spreading force upon droplet impact for the 
sprays’ stagnation flow field. In standard (0=0° orientation) 
spray cooling configurations the value corresponds with the 
heat flux (i.e. as the We increases, the heat flux increases to a 
critical value). The spreading force is felt normal to the heat 
exchange surface. In the case of angular trajectory sprays, the 
velocity vector must be decomposed whereby the vertical 
component (normal to the heater surface) is used in the We 
calculation. 


We = PyJ&L 

a 


(4) 


Table 3 gives a listing of the dispersion velocity and We for 
each 0 angle tested. The vertical velocity components decreased 
with increasing 0 as expected. The general trending 
correspondence for the We is inverted since the 0=30° and 45° 
orientations outperformed the 0° orientation (with exception to 
the straight fin "p0°, 0=45° case). 

One feature of the center weighted nozzle array used was 
that it created excess liquid in the center of the heat exchange 
surface for the 0=0° orientation. This was witnessed for both the 
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flat and straight fin surfaces. The straight fin structures 
themselves exacerbated the pooling and helped to make the 
excess liquid line more pronounced. As heat loads were 


e 

Vertical 

Velocity 

(m/s) 

Horizontal 

Velocity 

(m/s) 

We 

0° 

21.2 

0.0 

1796 

o 

O 

18.4 

10.6 

1348 

45° 

15.0 

15.0 

898 


Table 3 Dispersion Velocity and We for variable 0 values 

gradually increased towards CHF, the liquid line receded 
slowly. In the flat surface case, it was only slightly noticeable at 
CHF. However, in the case of the straight fins, it never fully 
receded. For the angular trajectory orientations (0=30° and 45°), 
the excess liquid line was not noticeable throughout any of the 
heat flux regimes. This along with the increased pressure 
coefficient indicates that the angular trajectories help promote 
better fluid spreading across the heater surface and ultimately 
better heat transfer. This is instituted by the sweeping motion of 
the spray droplets’ horizontal velocity components. Based on 
the experimental data, it appears that instituting a sweeping 
motion, while maintaining a dominant vertical velocity 
component fosters higher heat fluxes than the case of a 
dominant sweeping velocity. This is witnessed in the 
comparison of the results for the 0=30° and 45° orientations for 
each case. However, further investigation is required. 

The straight fin surface exposed to angular trajectories 
creates two unique flow scenarios with special considerations. 
When examining fluid motion upon the surface, both the y=0° 
and 7=90° (broad-crested weir) case may be consider open 
channel flows. For the 0=0° orientation testing, fluid was 
witnessed jetting out of the structure channels. This channel 
jetting motion was also witnessed in the y=0° case (applied to 
both 0=30° and 45°). However, it was negligible in the 7=90° 
case for both 0 =30° and 45°. Upon droplet impact on the 
surface, the non-evaporated bulk fluid must conform to the 
surface structures and the channel paths in its motion. This is 
inertia driven. As such, excess liquid in the 7=90° case 
undergoes deflection by the structures and redirection down the 
channel. This translates into lower channel velocities in relation 
to the 7=0° case. The lower channel velocities allow for fuller 
nucleation and ultimately higher heat fluxes. The higher channel 
velocities inherent to the 7 = 0 ° cases suppresses nucleation more 
and lowers the heat flux. 

SUMMARY AND CONCLUSIONS 

Spray cooling heat flux measurements were performed on a 
straight fin enhanced surface and a flat surface using PF-5060. 
Tests were performed under nominally degassed (41.4 kPa) 
environmental conditions. The volumetric flux (0.016 m 3 /m 2 s) 
and nozzle height from the surface (17 mm) were held constant 
for all the tests. 


The maximum heat flux attained occurred with the straight 
fin surface (Is) at 7=90°, 0=30°. The explicit value for the heat 
flux was 140 W/cm 2 giving an enhancement of 75% relative to 
the flat surface case at the vertical spray orientation (0K)°). 

The We decreased with increasing angular trajectory (0). 
Both angular trajectories of 0 =30° and 45° attained higher CHF 
values than the 0=0° case. 

The maximum heat flux occurred at a spray trajectory of 
0=30° for each of the case studies performed. This angle 
presented a dominant vertical velocity component in 
comparison to the sweeping (horizontal) velocity component. 

Application of the angular trajectories eliminated liquid 
pooling on both the straight fin and flat surfaces. 

The 7=0° case for the straight fin angular trajectories held 
higher channel velocities than the 7=90°. This corresponded to a 
higher level of nucleate suppression and lower heat fluxes. 
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